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Several methods of solution making use of a beam, plate or shell theory 
are available for the solution of straloht rntatlna. aswmmetrlr cross section 
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blades (refs. 1 to 3). The coupled equations of motion of such blades based 
on either linear theory (refs. 4 to 6) , or geometric nonlinear theory allowing 
for small precone (refs. 7 to 9, to mention a few) are available. However, the 
equations of motion Including large variable sweep for blades of advanced 
turboprop type configurations are not yet available. While finite-element 
modeling of the turboprop blades appears to be the most appropriate method for 
blades of such complex geometry, such studies with the existing codes at the 
Lewis Research Center revealed that the predicted results are satisfactory only 
for the first few modes. Furthermore, the complicating effects Included In the 
finite element codes, and also In the plate and shell theories, make the under- 
standing of the Individual and collective effects of the governing parameters 
Impossible. In order to conduct parametric studies to assess the various com- 
plicating effects, and to acquire a physical understanding of the complex blade 
dynamic problem. It Is proposed to use a simpler beam theory to model the 
rotating blade with the complicating effects successively taken Into account 
to reveal the relative Importance of the Individual and collective effects. A 
preliminary study made by using a set of linear equations of motion of a 
torslonally rigid, pretwisted, rotating blade Including Coriolis effects was 
reported In reference 10 wherein the effects of sweep on the dynamic behavior 
were Introduced by preconing the blade with respect to the plane of rotation. 
The effects of linear pretwist, precone and linear Coriolis effects the 
vibration and stability of rotating blades were discussed In reference 10, and 
It was pointed out that the Coriolis effects must be Included In the analysis 
of thick blades but could be disregarded In analyzing thin blades possessing 
small pretwists. The position and width of an Instability region was shown to 
be dependent on the extent of pretwist, precone and whether or no l the Coriolis 
effects were Included In the analysis. Although the blade was considered to 
be torslonally rigid, and hence the results somewhat restricted In generality, 
considerable Information on the various governing parameters was obtained. The 
first objective of the present effort Is to determine the effect of second 
degree geometric nonllnearltles on the natural frequencies, steady state 
deflections and mode shapes of the blade cases considered In the previous 
Investigation (ref. 10). The second objective Is to find the parameter limits 
within which the second degree geometric nonllnearltles are adequate to 
properly represent the blade dynamic characteristics, by comparison of results 
produced by beam theory to those produced by MSC NASTRAN. It may be noted here 
that only the second degree geometric nonlinear effects are Included In the 
present beam theory together with Coriolis effects. Further, there Is no 
restriction on the degree of nonlinearity In the MSC NASTRAN although the 
Coriolis effects are not accounted by this finite element code. Thus, a fair 
comparison of frequencies and steady state deflections produced by the present 
beam theory to those from MSC NASTRAN (for such blade configurations that are 
Insensitive to both linear and nonlinear Coriolis forces) would establish the 
validity of the restriction of the nonllnearltles to only the the second 
degree. Further complexities of torsional, extenslonal, rotational and warp- 
ing couplings can then be addressed once the accuracies of the present 
restricted beam model are properly validated. 

In order to accomplish the stated objectives, the required equations of 
motion are derived by using the theory presented In reference 9, and by retain- 
ing geometric nonllnearltles up to second degree. The Galerkln method, with 
nonrotating normal modes. Is employed for the solution of both steady state 
nonlinear equations and linearized perturbation equations. Parametric 
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EQUATIONS OF MOTION ANU METHOO OF SOLUTION 

The coupled f lap-lag-extenslonal equations of motion of a rotating, tor- 
slonally rigid, linearly pretwisted and preconed blade of uniform rectangular 
cross section, shown In figure 1, Including Coriolis effects and second degree 
geometric nonllnearltles but disregarding all other higher order effects, can 
be derived by using the theory presented in references 9 and 11. Such equa- 
tions are presented below (a list of notation Is given In appendix B): 

Flatwise bending: 

mw - mQ 2 w sin 2 B pr + 2mQv sin B pc ♦ ntt 2 sin B pc cos B pc fx + u - u f ) - { T w* ) ' 

+ w" (El cos^ e + El sin 2 e) ♦ v B (EI - El ) sin e cos e> = 

I r\r\ 55 nn ) 

\ 

- jft 2 sin B pc cos B pc p (1^ sin 2 e ♦ cos 2 6)| (1) 

Edgewise bending: 

mv - m£l 2 v + 2mQ cos B pc (u - Up) - 2mOw sin B pc - (Tv 1 ) 1 

+ < w B ( EI__ - El ) sin e cos e + v"(EI sin 2 e + EI PP cos 2 e» = 

I ' EE nn nn EE ) 

- in 2 sin B_ cos B„ p (l ct - I ) sin e cos e( (2) 

) Pc Pc EE nn \ 

Extension: 

m[ u - u F - ft 2 cos 2 B pc (u - Up + x) + wQ 2 sin B pc cos B pc 

- 2v« cos B pc ] - ( AEu 1 ) ’ - 0 (3) 

where 

T - - ^ m[(i - \ip - G 2 (R + x - Up) cos 2 B pc - 2M cos B pc 


♦ fl w cos B pc sin B pc 


]dx, 


( 4 ) 
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u. = \ f* (V 2 + w' 2 )dx, 
r 5 


( 5 ) 


and 


m 


//pdy dz, » =// dy dz, l 5? *// y2 <•» dz - : nn \f/ z2 dy dz ' 


( )' " fc ()•<’)■ it < ) (6) 


Defining the following parameters, 

w = w/L, v = v/L, i) = x/L, t = SJt , R = R/L, etc., 


(7) 


assuming solutions are separable In time and space, and making note of the 
following relations 

T- = 7-*T l = fs'*S? = Q H" etc - (8) 

dx dn dx L dn dt Qt 

one can rewrite equations (1) to (3) In the following nondlmenslonal forms: 

* l n p 2 

w ♦ 2 sin B pc ^ - W sin 2 B p{; - \ sin 0 pc cos B pc f (v 1 ♦ w' )dn 


- cos' 0 pc (w" Q - w' S) - 2 cos B 


pci*" J 5 

L n 


v dn - w'v 


sin B pc cos B pc |w" J w dn - w'wj + *w iv ^ 


p h 2 

cos' e + ^ sin e 
a ) 


r. 1V,1 


+ w h, (2y5 Sin 26)1 - 1 I + w"(2y's cos 26) ^ - 1 ♦ ?v (g stn 26)1 - 1 


+ v"'(2yS cos 26)^ - lj - v"(2 Y 2 5 sin 26)^ - 1 J * s1n B Pc cos B Pc Q 

1 1 p 1 2 

+ w" | a dn - W" j ii cos‘ B pc dn - w'u + w u cos 0 pc 

n n 

- -n sin B pc cos 0 Rc - ^ s1n B Pc cos B Pc s1n 20 


(9) 
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V - 2 sin 0 Dr w - v - cos 2 0 n „(v"Q - v'S) - 2 cos 0„.|"v" f v dn - v'vl 


rL L J 


- cos 


B Pc d T j" ( v ' 2 + w ' 2 ) dri + s1n B p c cos B p c £ v “ J W dn - v'wj 


+ w 1 ^ sin 2ej^ - + w" 1 ( 2y£ cos 26)/^ - 1^ 


) 


/ 2 \ ( 2\ / 2 

- w"(2t 2 5 sin 20)^ - lj + v 1v ?|s1n 2 6 + ^ cos 2 A - v"'(2y 5 sin 2e)|^ - 

- v"(2y 2 5 cos 2Q)(^ -1^ + 2 cos 0 pc u - v" j” u cos 2 0 pc dn 


I Y 


+ v" | u dn - v'u + v ! u cos 2 0 pc = - sin 0 pc cos 0 Rc ( cos 26 

A L 




(V' 2 + w' 2 ) d n 


- cos 2 0 n J + i cos 2 0 pc f ( v ' 2 + w' 2 ) dn 


( 10 ) 


Pc 2 


+ w sin 0 Rc cos 0 pc - 2v cos 0 pc - (AE/mQ 2 L 2 )u" = n cos 2 0 pc 


(11) 


where 


0 » R(1 - n) + 0 


•5^1 - n 2 ). 


S = (K + n), (w)' = S_ (w), (w) = t- (w), 


dT 


5 = (EI^/pAl/fi 2 ), and 0 = <p r Y n (12) 


Before discussing the method of solution, It Is worthwhile to point out 
the various Important linear terms associated with precone, and also the non- 
linear terms existing In the present equations. The linear terms associated 
with precone are addressed first. Referring to equation (9), one can see that 
a linear softening term, (-w sin 2 0 Pc ), appears In the flap equation. This 
term vanishes for zero precone, but becomes an Important term for suitably 
large values of precone and rotational speeds, and contributes to the mechanism 
of rotation Induced Instability. Next, the terms (2 sin 0 Pc 7) and (-2 sin 
&Pc 1" equations (9) end (10) respectively are the linear Coriolis force 
terms that arise due to the Inclusion of precone. The effect of these terms 
on the linear frequencies has been discussed In detail In reference 10. Con- 
sidering equation (11), one observes that there Is one linear term (w sin 0p c 
cos 0 Pc )jvh1ch vanishes for zero precone, and that the linear Coriolis force 
term (2 v cos 0p c ) appears In this equation whether or not precone Is present 
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An +he ^privation of the equations. Since It will be shown that the Inclusion 
of extenslonal degree of freedom Is not that Important tor analyzing ooui tn... 
and thick blades, further discussion of the other nonlinear terms assoc lated 
with the extenslonal deformation will not be attempted In this section. Next 
the Important nonlinear terms existing In equations (9) and (10) are considered. 
The nonlinear Coriolis force term's are shown by underlining the m once, and 
these terms will be present In the equations even when precone Is ? bsent ; Tb * 
nonlinear terms which are shown by underscoring them twice In equal Ions (9) a d 
(1G) are the contributions from the tension terms (Tw ) and ( Tv ) • ^s are 
of the three nonlinear Coriolis terms just discussed. Howeve. , the d0 “ bl J 
underlined terms vanish for zero precone. Finally, the ter |" shown by bro 
underlining Is the effect of foreshortening of the blade This term also va - 
ishes for zero precone. It may also be mentioned here that for the 1 siting 
case value of 90° precone, all nonlinear terms In equations (9) to (11) vanish 
excep^nS those associated with extenslonal Inertia (u, u F . Since the non- 
■Mnear terms associated with extenslonal deformation are already noted to be 
the effect of geometric nohlloearltle, should become almost neo- 

llglble at 0p c = 90°. 

ThP f laD-laa-extenslonal equations are solved by the Galerkln method by 
assuming that the dimensionless bending and extenslonal de f n erms 

of a series of generalized coordinates and mode shape functions are as fo 


where 


w 

II 

,w o3 * 

AWj (T))^j(n) 

(13) 

V 

II 

-M 

(v oj * 

AVjM^n) 

(14) 

u 

II 

(u oj * 

AUj(x))ej(n) 

(15) 

(6 3 n) 

- COS 

(Bjn) 

- a^[s1nh (Bjo) - sin (B^n)] 

(16) 


0j(n) = 2 

sin (T^n) 

(17) 


y 

3 •'« 

-b 

(18) 


Equations (16) to (18) are the nonrotating normal modes for a cantilevered beam 
n“d at ; . 0, and free at , .1. furthermore, the quantities «o1. *03 «"<• 

U-. in the generalized coordinates constitute the equilibrium quantities 
while aw j , Avj and Auj are the perturbation quantities. 

By substituting only the steady state equilibrium quantities Into * h ® 
nonlinear equations (9) to (11), assuming n-normal modes each of the vari- 
ables u, v, and w and carrying out the Galerkln process tradltlona ■ ly , one 
obtains a set of 3n nonlinear equations In terms of w oj* v oJ a " d u oJ- The 
constants a\ and Ba are taken from reference 12, and the resulting 
equilibrium equations are solved by using a computer program based upon a 
finite-difference levenberg-Marquardt algorithm (ref. 13). Next, 
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nm.a+inn* m> to (15) are substituted into equations ( 9 ) to (n), the Jalerkln 

msm&S^Maotr^ 

matrix notations: 


[L]t X o> + CNLllXo 1 - (S) 


[H]{X) + [C]{X} + 1*1 W ■ 0 


(19) 


( 20 ) 


where 


* (w ol , w o2 w on , v or v o2 , .... v on , u ol . u o2 


u )T 

on 


( 21 ) 


X = (aw r aw 2 , 


AW n , AV r AV 2 , 


AV p , AU r 


AU 


?’ 


«,.iT 


( 22 ) 


with t and LHbelngrespectWelythe lj"®ar and an q stlf f- 

rotating blades. 


RESULTS AND DISCUSSION 

The nonlinear steady state -!. 1 ( 2 ; ) a wJ r ^!olve^by 1 u^ng°cSuter t 

KJTsSKTp/SJ -Tn^e presence of '-^^^""r^^er^irs^e™. 
will occur In pairs of P“J e ']L'™S ?4quency equation (20) reduces to a 

general computer program Typical th'ckn. s rat » „ (d/b , 0.25) 

sent advanced turboprop blades (d^ - u.u^, or p v 1Q Jt ^ 5e noted 

that the* radius or^dls^Y Is assumed to be zero for simplicity. 
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Convergence 


The convergence of solutions produced by the Galerkln method for the 
coupled flap-lag-extenslon equations, using various numbers of nonrotating 
normal modes for the Independent variables, Is shown In table I. The blade 
considered for this convergence study has a precone of 30° and a thickness 
ratio of 0.5. The blade chord at the root Is set perpendicular to the axis of 
rotation (<p * 0°) and the blade rotational speed Is one half of the funda- 
mental mode frequency of the same nonrotating blade (Q/w-| = 0.7/. Natural 
frequencies for this blade with zero pretwist are determined by varying the 
number of nonrotating normal modes, n, used In the series solution assumed. 

It can be seen from this table that a five-mode solution produces the frequen- 
cies of the linear equations, and also those from a perturbation solution of 
the nonlinear equations, converged tc five significant figures. A further 
comparison of the present flap-lag-extenslon frequencies, (p/\-|), obtained 
from the solution of the linear equations to those giver In references 10 and 
14 shows an excellent agreement. Further results for various combinations of 
pretwist, precone, setting angle and rotational speed are obtained by using a 
five-mode Galerkln solution, and the Individual and collective effects of the 
various parameters are discussed In the following sections. The validity of 
restricting the geometric nonllnearltles to second degree only Is assessed by 
comparison of the present beam theory results to those produced by MSC NASTRAN 
for specialized cases of thin blades In the following section. 


Comparison of Present Results 

Comparison of the frequencies from the solution of the present linear 
equations and those from a perturbation solution of the nonlinear equations Is 
made to the frequencies produced by MSC NASTRAN using 250 CQUAD4 elements In 
tables II to IV tor typical values of precone, rotational speed and thickness 
ratio. The steady state equilibrium deflections produced by the present beam 
theory and those produced by MSC NASTRAN are compared In figure 2. Considering 
the results presented In table II corresponding to a thin blade poaching zero 
pretwist, zero precone and zero setting angle, one can see that the lowest six- 
mode frequencies obtained from perturbation solution of the beam theory equa- 
tions agree to within one half of 1 percent with those given by MSC NASTRAN, 
for wide range of rotational speeds. It may be noted that the present flap- 
lag-extenslonal equations cannot predict the torsional frequencies since this 
degree of freedom Is not considered In this study. Next, a comparison of the 
lowest three frequencies of the same blade considered earlier but with a 15° 
precone are presented In table III. A further comparison of linear and non- 
linear frequencies Is also made In this table. Here also, the agreement 
between the two sets of results Is good. The effect of geometric nonllnear- 
ltles on the lowest three flatwise modes Is seen to be of a stiffening char- 
acter, and the frequencies are found to Increase with Increasing rotational 
speeds due to the nonllnearltles. Table IV shows further comparison of results 
for large values of precone, for both thin and thick blade cases, at various 
rotational speeds. Considering the trend of results observed so far for thin 
blades, It Is evident that the first bending modes In both flatwise and edge- 
wise directions are affected much more than higher modes. The percent dif- 
ference between beam theory results and MSC NASTRAN results Increases with 
Increasing precone for a given rotational speed, and for a given precone with 
Increasing rotational speeds. However, for practical rotor speeds of the order 
of (fl/w|) <1, the difference between the two sets of results Is not 
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greater than 5 percent. When a comparison of the results obtained for thick 
blades (d/b = 0.25), Is made, It appears that the percent error observed for 
the second mode Is quite large. In order to assess the effect of Ignoring the 
Coriolis effects, results for the thick blade cases are also obtained from the 
present beam theory by Ignoring the Coriolis effects, and these results are 
compared to the corresponding ones from MSC NASTRAN. These are also Included 
In table IV. From thl. comparison of results. It Is evident that the results 
from MSC NASTRAN are In closer agreement with the corresponding results from 
the beam theory when the Coriolis effects are Ignored In beam theory. It 
should be noted here that the rather large difference observed for the second 
mode frequency of thick blades, (which Is the fundamental edgewise frequency), 
may be partially attributed to the fact that the present slender beam approxi- 
mation In the beam theory may be Inadequate to predict the frequencies of 
stubby blades In general and the stiff edgewise mode frequencies In particular. 
Proper comparisons of these frequencies Is only possible by Including the 
shear and rotary Inertia effects In both beam theory and MSC NASTRAN. 

Next, the steady state deflections for thin blades given by the present 
beam theory are considered. For the untwisted blade cases considered In this 
work for various precones and rotational speeds. It Is observed that the flat- 
wise deflection Is the most significant while the edgewise and extenslonal 
deflections were almost Insignificant. The extenslonal deformation was found 
to be more significant than the edgewise deflection. However, for pretwisted 
blades, both flatwise and edgewise steady state deflections became quite sig- 
nificant, and the magnitudes of these deflections were found to be large In 
comparison to the corresponding untwisted blade deflections. Thus, for the 
untwisted, thin blade cases, the distribution of dimensionless flatwise 
deflection, w, along the length of the blade obtained from the present beam 
theory Is compared to the corresponding one from MSC NASTRAN In figures 2(a) 
to (d), for typical precones and rotational speeds. An examination of these 
results Indicates that the trends shown by both sets of results considered are 
consistent. Furthermore, the agreement between the two sets of results Is 
extremely close for low rotational speeds. The difference between beam theory 
results and MSC NASTRAN gradually Increased from zero at the root section to 
the maximum at the tip In general. For precones of the order of 45° and 
rotational speed parameter value of up to (8/u-|) = 1.0, the great st difference 
between the two sets of results, (steady state deflections and frequencies up 
to third mode), has been found to be of the order of 5 percent for thin blades. 
It Is Interesting to note that the deflections given by the beam theory are 
consistently greater than those produced by MSC NASTRAN. 

From the foregoing discussion of results, It appears that the present beam 
theory, Including geometric nonl Inearltles up to second degree only, predicts 
the natural frequencies and steady state deflections to an acceptable degree 
of accuracy In the case of thin blades having precones of up to 45° and for 
blade rotational speeds that are practically encountered In their applications. 
It Is believed that the Inclusion of the torsional degree of freedom, which Is 
necessary for calculation of stability boundaries for the blades, should not 
alter this trend of agreement of results, since the torsional mode coupling may 
affect the stability boundary but not the convergence of the first few modes 
that are well separated from the basic torsional mode frequency. It may thus 
be concluded that the second degree nonl Inearltles are adequate for modeling 
thin blades with large precone and for rotational speeds encountered In their 
practical operational range. 
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In order to understand the Individual and combined effects of ■ 
linear and nonlinear Coriolis forces, various nonlinear terms (seco I .fgree 
geometric nonllnearltles) and pretwist on the frequencies of rotating 
parametric studies were conducted for two typical thickness ratios of c . 

0.25 for various setting angles (collective pitch). These results aro pr*- 
sented In tables V to IX. Figure 3 shows the effect of setting angle, pretwist 
and precone variations on the fundamental mode frequency parameter, (p-j/X-j). 


Effect of Varying Precone 

In order to determine the effect of varying precone on the frequency par- 
ameter ratios of untwisted blades with zero setting angle (<p = 0), the flap- 
lag-extenslon equations were solved for a typical rotational parameter value 
of (Q/w-]) = 1.0. The value of precone was changed from 10° to 50°, and the 
linear and nonlinear frequencies were determined which Include the Coriolis 
effect terms also. These frequencies are listed In table V together with the 
percent frequency variation based upon the nonlinear frequency for each mode. 
The following observations are made from the results presented In table V. 

1. For a given rotational speed and thickness ratio, the flatwise mode 
frequencies decrease with Increasing precone In the case of both linear and 
nonlinear theories. The first edgewise mode frequency, (refer the second mode 
frequency of thick blade with d/b = 0.25), given by linear equations or non- 
linear equations shows an Increasing trend with Increasing precone. 

2. The flatwise mode frequencies produced from the perturbation solution 
of the nonlinear equations are higher than the corresponding frequencies 
obtained from the linear set of equations. The fundamental edgewise frequency 
given by the solution of nonlinear equations’ Is lower than the corresponding 
linear solution value. 

3. The effect of geometric nonllnearltles, as could be seen from the per- 
cent frequency variation, Increases with Increasing precones of up to 50° con- 
sidered here. The frequency change for flatwise modes Is seen to be positive 
(stiffening) while for edgewise modes it Is negative (softening). However, It 
should be noted that for 90° precone, the effect of geometric nonllnearltles 
becomes almost zero as should be expected. 

4 . The fundamental mode frequency shows the strongest frequency variation 
(positive for flatwise mode and negative for edgewise mode) due to the presence 
of nonllnearltles, and this effect Is se^n to decrease as the mode number Is 
Increased. 


Effects of Pretwist, Coriolis Forces and Geometric Nonllnearltles 

The Individual and combined effects of pretwist, precone, Coriolis forces 
and thickness ratio on the frequencies of rotating blades In the absence of 
geometric nonllnearltles were presented and discussed In reference 10. How- 
ever, the linear frequencies are presented again In all the following *ables 
for the purpose of completeness, and also to provide an easy access fo 
discussing all the results together. The slight differences (at the lourth 
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;4nn<fAr 2 nf f^nnro) one may find between the linear frequencies presented In 
this paper and^hose In reference 10 are due to the fact that the extenslonal 
degree of freedom was Ignored In reference 10 while this degree of freedom Is 
Included In this work, and also due to the fact that the methods of solution 
used In these two works are different. 

In order to ascertain the Individual and combined effects of precone, 
rotational speed, Coriolis forces and the various terms that arise In the 
equat.ons due to the Inclusion of second degree geometric nonllnearltles, 
several cases rotating blades were solved. These results are presented In 
tables VI and VII for both untwisted and 30° pretwisted blades having zero 
collective pitch (<»> * 0®). Typical rotational parameter values of (G/o>i) - 0.5, 

0.8 and 1.0 were considered In this study since these values generally encom- 
pass the practical operational speeds of advanced turboprop blades. The fre- 
quency parameter ratios (p/X]) obtained from the solution of the linear 
equations, Including or excluding the Coriolis effects, are presented first In 
these tables. Next, the frequencies obtained from the perturbation solution 
of the nonlinear equations are presented, starting with the frequencies of full 
nonlinear equations followed by those obtained by Ignoring one key parameter 
In the nonlinear equations at one time. Tnus, the frequency parameters under 
the column with = 0 represent the frequencies obtained by Ignoring the 

foreshortening effects although all other effects are present, those under the 
column with (D^|< * 0, = 0) Illustrate the effect of Ignoring the nonlinear 

terms arising from the centrifugal tension terms ( ( Tw ' ) and (Tv') 1 ), while the 
frequencies In the last column Indicate the effect of Ignoring the nonlinear 
terms associated with the extenslonal deformation together with these from 
centrifugal tensions. While the results presented In this form are useful for 
future comparison, the Individual effects will be clearer If percent variation 
of the frequencies are calculated for each category based on full nonlinear 
solutions. This Is accomplished In tables VIII and IX. The following observa- 
tions can be made from the results presented In these tables: 

1. The limiting values of precone of zero and 90° are considered first on 
the vibrational characteristics. When Bp c = 3, one can see from equations 
(9) to (10) that the flap-lag equations are coupled through the nonlinear 
Coriolis force terms together with the linear and nonlinear extenslonal 
deformation coupling terms. Coupling due to the latter category of terms Is 
not that Important however. Pretwist In the blade brings In the additional 
Important structural coupling between flap and lag deflections. Thus, for 
untwisted and pretwl'ted blades, the effects of geometric nonllnearltles 
(excepting for the Coriolis terms which are Important for thick blades) on the 
natural frequencies Is almost negligible when precone Is zero. When the pre- 
cone Is 90°, the flap and the lag equations are coupled through linear Coriolis 
force terms and extenslonal Inertia even for the untwisted case, while the 
extenslonal equation of motion Is coupled through the Inertia associated with 
foreshortening. Since the right hand side for this case of 90° precone Is zero 
for all the coupled equations, the steady state deformations will be absent, 
and the equations produce results that depend only on Coriolis effects. This 
can be verified from the results presented In table VI(a) for 90® precone case. 
The effects of geometric nonllnearltles and linear and nonlinear Coriolis 
forces are therefore Important for precone angles other than these extremes. 

2. Although the frequencies show both Increasing and decreasing trends for 
a given precone with an Increase In rotational speed, (refer to table VI(a)), 
the effect of geometric nonllnearltles Is seen to Increase the frequencies of 
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flatwise modes, and to decrease the frequency of first edgewise mode (refer 

f ahl ftf WTTT a nrl T V \ f AIT r f aH a a H A r A 4-1 il r 4- nrl K1 f a r a r a a 4- A A a'i v>a^ a r 

www • v j • * * i vmu am/ i vi wi ii v w * ^ i, vvi wnvj ^i\.v.nijtv-u w I uu v, J I ui u I I p i uw I, I v.u I I u but 

speeds considered 

3. The effect of second degree geometric nonllnearltles Is the greatest 
on the fundamental mode, and decreases as the mode number Is Increased. The 
effect of second degree geometric nonllnearltles Increases d, the rotational 
speed Is Increased. By referring to tables VI I ( a ) and (b). It can also be seen 
that an Increase In precone from 15° to 45° Increases the effect of the geo- 
metric nonlinear terms considerably for a given rotational rpeed. 

4. In tables VIII and IX, the percent frequency variation due to the 
absence of linear or nonlinear Coriolis effect terms In the presence the other 
and with second degree geometric nonllnearltles Is presented. These frequency 
variations are calculated based upon full nonlinear frequencies. It can be 
seen from these results that the nonlinear Coriolis effects are stronger than 
the linear ones In affecting the frequencies of all modes for high rotational 
speeds generally. Both hardening and softening characteristics are exhibited 
by the linear and nonlinear Coriolis terms for a given mode with Increasing 
rotational speed. This trend Is to be attributed to the fact that even though 
the linear (or nonlinear) Coriolis terms are Ignored and nonlinear (or linear) 
Coriolis terms are retained In the equations, the presence of other geometric 
nonllnearltles also affect the resulting frequencies which causes these mixed 
trends. It emerges clearly, however, that both linear and nonlinear Coriolis 
effects are Important for high thickness ratio untwisted or pretwisted blades, 
and these effects are Insignificant for low thickness ratio blades. For the 
cases where the Coriolis effects are Important, the linear Coriolis effects 
oppose the nonlinear Coriolis effects. The Influence of the linear and non- 
linear Coriolis effects are the greatest on the first flatwise mode and on the 
first edgewise mode, and are Insignificant on other higher modes. These 
effects are more pronounced for larger precones and higher rotational speeds. 
One may thus conclude that both linear and nonlinear Coriolis effects can be 
Ignored In analyzing thin blades which are typical of advanced turboprop blade 
configurations . 

5. The effects of Ignoring either the foreshortening terms, or the 
centrifugal tension coupling terms on the nonlinear frequencies are shown In 
the last two columns of tables VI J I and IX. It can be seen from these results 
that these two effects produce the greatest variations on the frequencies, and 
that the first mode Is affected to the greatest extent. It may also be noted 
that these terms arise due to the presence of prerone, and the frequency 
Increases by nearly 20 percent due to the presence of the second degree geo- 
metric nonlinear terms. 

6. Tne onset of static Instability for various cases of preconed rotating 
blades with and without pretwist was predicted by using the linear equations 

In reference 10. It was shown that a 60° preconed blade with a thickness ratio 
of 0.05 become., statically unstable for 1.48 < (Q/o>-|) < 1.49 If linear 
equations Including Coriolis effects were used for the prediction of the 
Instability. By using the present second degree geometric nonlinear equations 
with Coriolis effects, this Instability was found to occur for 1.13 < (ft/u-|) < 
1.14 for an untwisted blade, and for a 30° pretwisted blade. When the 
untwisted thin blade with 60° precone was solved by using MSC NASTRAN, It was 
observed that the pseudo static configuration became unstable at (Q.'w|) = 0.8. 
Since the results from the analysis using MSC NASTRAN gave consistently good 
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agreement of frequencies 

ones from the present Mi ^^l^n^rcoupil^ (which Is present In 
linearities, It Is believed that theorv) mus+ have been responsible for 

HSC NASTRAN analysis bu t absent 1 bea, ■ « that the torstonal 

S^p’llng'ltt MS* m the beam theor, for a fair prediction of 

Instability boundaries. 

7. The effect of pretwist -“P’^^'rtslSns^l^r^lso 
blades was studied In refeience 10. Sli c nc non ^oear Coriolis 

since the second degree geometric nonl Inearl ties ana xne the 

effects have not affected the higher modes ^Kd cSSpl^ng trend of 

coupling trends for thin blades. The w ^ wU * a mode frequen cy) and Increasing 
decreasing the lowei frequency (firs g . f the tw0 c i 0 ser modes 

the higher frequency (second flats, se ^ e "n the results presented In 

of untwisted blade due to ’ * v !"e for the precones and rota- 

tables VlII(a) and (b) for the thick blad , geOTle tr1c nonllnearlt- 

»‘K?iltr?r« « the frequenc1es S of°pretw1sted^blades are similar to 

those observed for untwisted blade cases. 

8. frequency parameter ratios found^hatlorT;! 

from the present linear and nonlinear equ • • nonlinear frequen- 

blade (d/b = 0.05), having 0 Pc = « and rS* the respective linear f re- 
cles of the lowest three modes ^re grea These frequency variations 

quencles by about 0.067, 0.005 and 0.00 J of ? 979 , _ 0 .263 and 0.137 

for a thick blade (d/b - °-25) w these results with those obtained for the 
percent respectively. By comp Scented In table VIII(b), one can conclude 
zero setting anlge case (j - I n^rU^ are far more severe for «, - 0* 
that the effects of geometric . «qo f or nretwlsted blades 

than for * - 90V Results pertaining^ * - «0 for P retwi s dUcussed 

also show similar trends to those observed J" ™J d “ n % lnce U has been 
above. F or brevity these resu nea r1t1e S and Coriolis forces affect the 

established that the geometric nonl 1" ^Undesirable to present the variation 
fundamental mode most severely , ^ ^ . t1o wUh respect to the variation 

of the fundamental mode frequency param eter ratio wv cn re This Is 

of precone for various combinations 0 th « s f^ qur e that for a given rotational 
shown in figure 3. It can be seen from th s flflur that mode fre . 

speed and precone, the variation of ^etwmchang s ^ 458 

quency to an appreciable exten distinct and well separated for each 

the fundamental mode frequency Is quite distinct anowei p ^ 

combination of setting ««’V ™1 ^ ^ SreSne. Since for Bp c - 90-. 

quency values droop down to a s ch a ft the effect of setting angle 

tt- nreconed blade b ®^ e s a cantlle^^d^h ft.^ frequency values must be 

vanlsno, and the slight dlf ar uina from pretwist and the linear 

attributed to the coupling obtained for the various cases discussed 

Coriolis terms. Finally, results ^ , . . The differences observed were 

earlier by Ignoring tbe compared to the corresponding 

In the fourth or fifth For geometric and physical par- 

flap-lag-extenslonal Motions 5 ?]“ e xtenslonal degree of freedom can thus be 

ameters considered in ths work t J e shapc s calculated about the 

safely Ignored. It Is also ODserveo ^ the present nonlinear 

tpreclable extent from those obtained from the 

linear theory. 
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CONCLUDING REMARKS 

Coupled f lap-lag-extenslonal equations of motion of rotating, pretwisted 
cantilever blades of uniform rectangular cross section are derived Including 
large precone, Coriolis effects and second degree geometric nonllnearltles. 
Parametric studies are conducted to assess the Individual and combined 
Influence of the various complicating effects by solving the nonlinear 
equations using a linear pertubatlon technique. The following major 
conclusions have emerged from the present effort: 

1. Inclusion of the geometric nonllnearltles up through second degree 

only In the flap-lag-extension equations appears to be adequate for the pre- 
diction of steady-state deflections ana the first few natural frequencies (If 
they are well separated from a basic torsional frequency) of thin blades hav- 
ing precones of up to 45°, and rotating at speeds of the order of = 1.0. 

The second degree geometric nonllnearltles show a stlfflnlng effect on the 
flatwise modes In general, and a softening effect on the first edgewise mode, 
for precones of up to 50° and for all rotational speeds considered In this 
work. The effect of geometric nonllnearltles on higher modes Is not signifi- 
cant. The greatest effect Is from those nonlinear terms which vanish for zero 
precone. However, If precone Is substantial, the second degree nonlinear terms 
can produce frequency changes of engineering significance. The Increase In 
frequencies of first flatwise and first edgewise modes are typically of the 
respective orders of +20 percent and -4 percent for blades with 45° preccne and 
for a rotational speed parameter of fl/u-| = 1.0. 

2. The effect of nonlinear Coriolis forces Is most severe on the first 
flatwise and first edgewise modes, and Insignificant on higher modes. The 
effect of Coriolis forces Is found to be significant for thick blade cases. 

In general, the nonlinear Coriolis forces oppose the linear ones, the nonlinear 
effect being stronger. Thus, both linear and nonlinear Coriolis effects can 
be Ignored In analyzing thin blades which are typical of advanced turboprop 
blade configurations. 

3. Preconing has significant Influence on the first flatwise and edgewise 
modes. An Increase In precone at a given rotational speed shows a softening 
effect on flatwise modes generally, and a stiffening effect on the first edge- 
wise mode. However, the softening effect Is much more pronounced than the 
stiffening effect. 

4. The coupling trends for pretwisted blades observed from the solution 
of the present nonlinear equations, do not differ to any appreciable extent 
from those ob c erved from the linear theory. This may be attributed partially 
to the fact that geometric nonllnearltles do not significantly affect higher 
modes for the present f lap-lag-extenslonal equations. 
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APPENDIX A: THE GALERKIN INTEGRALS AND MODEL EQUATIONS 

The various Integrals arising from the Galerkln process are defined 
below, and these are used In representing the modal equations In matrix forms 
subsequently: 


f' r 1 

l 1j = J Vj dn = j e 1 e j 


dn 


IJk 


/*’ 

| *J(x)t k (x) dx dn 


J « 


J dn 
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C 1J ’ 

1 dn 
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D 1Jk = j Vj S *k<*> d * dT > 
o n 

J M 

E uk = ) dn 
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■ j *>*] 


cos 2 0 + ^ s1n ^ d n 
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cos 26 dn 
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APPENDIX 8 - NOMENCLATURE 


A 1Jk* B 1J’ L 1» etc 

b, d 

d/b 


{B}, {X}, {X Q > 


cross-sectional area of blade 
modal Integrals (see appendix A) 

breadth (chord) and thickness of blade 

thickness ratio 

vectors 


[C] 

*** 

E 


I , I-- 


1, J, k 

L 

[L], [LN] 


m 

[M] 

n 


P 

R 

T 

t 


u, v, w 


modal damping matrix (gyroscopic matrix) 

Young's modulus 

area moments of Inertia about major and minor principal 
centroldal axes, respectively 

dummy Indices 

length of beam 

Linear and nonlinear components of the matrix representing 
steady state equilibrium equations 

mass of blade per unit length 

modal mass matrix 

number of nonrotating modes for each of the flap bending, 
lead-lag bending, and extenslonal deflections 

natural radian frequency 

radius of disc 

blaae tension 

time 

displacements of the elastic axis In X, Y, Z directions, 
respectively 


u, v, w 

V V W oj 

4 “r *Y 4w j 


dlmenslonaless deflections 
steady-state equilibrium deflections 

perturbation quantities 
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X 


y, * 

*r e r 

0p c 

Y 

*1J 

5 


n 

e 

M 

p 

♦ 

*j(») 


«l 


n 

( )' 
C) 


running coordinate along X-axis 

centroldal principal axes of beam cross section 

constants for assumed mode shapes 

precone angle 

total pretwist of the blade over Its length 
Kronecker delta 

nondlmenslonal rotational parameter, EI^/pAL^R^ 

nondlmenslonal length coordinate, x/L 

geometric pitch angle, <p + yn 

frequency parameter t> /EI nT ,/pAL* 

mass density of blade material 

setting angle (collective pitch) 

nonrotating flap and lead-lag bending mode shapes 

dimensionless time, Qt 

exact fundamental mode frequency of straight, nonrotating 
beam, 3.51602 \i 

rotor blade angular velocity, rad/sec 

primes denote differentiation with respect to x or n 

dot over a parameter represents differentiation with 
respect to t or t 
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TABLE I. - CONVERGENCE PATTERN OF FREQUENCY RATIOS (p /»i > 

OF A PRECONED, ROTATING BEAM PROOUCEO BY THE GAIERKIN 

METHOD WITH NONROTATING NORMAL MODES ORIGINAL PAfaF IS 

r«M . 0.5. 8 p c - 30*. d/b - 0.5, L/d . 20, ? - o*. OF POOR QUALITY 

Y - R - 0.] 


(a) Frequencies from solution of linear equations 
Coriolis effects 

Ignoring 

Mode 

n - 1 

n - 2 

n m 3 

n - 4 

n - 5 

1 

2 

3 

A 

3.788724 

7.008871 

108.883020 

3.788037 
7.008777 
22.355387 I 
44.204056 

3.787921 
7.008761 
12 . 355386 
44.204056 

3.787906 

7.008758 

22.355314 

44.204047 

3.787903 

7.008758 

22.355303 

44.204045 

5 


108.882980 ( 

52.025490 

62.025447 

62.025401 

(b) 

requencies fr 

om solution of linear equations in 
Coriolis effects 

eluding 

1 

2 

3 

3.6*5308 

7.302023 

108.92232 

3.634662 

7.301872 

22.331668 

44.247993 

3.634553 

7.301851 

22.331666 

44.247881 

3.634539 

7.301849 

22.331595 

44.247870 

3.634535 

7.301848 

22.331584 

44.247867 

5 

— 

108.92610 

62.017091 

62.017048 

62.017003 

(C 

) Frequencies from perturbation solution of nonl 
equations: Coriolis effects ignored 

Inear 

1 

2 

3 

A 

3.847417 

7.019543 

109.47317 

3.846369 

7.019438 

22.378138 

44.213271 

3.846218 

7.019422 

22.378096 

44.213273 

3.846199 

7.019420 

22.378012 

44.213263 

3.846194 

7.019419 

22.377999 

44.213262 

H 

b 


109.57995 

61.959160 

61.958438 

61.958267 

(< 

i) Frequencies 
equati 

from perturbation solution of non 
ons: Coriolis effects Included 

linear 

T“ 

2 

3 

A 

3.742021 

7.214380 

109.51676 

3,740674 

7.214485 

22.363212 

44.242350 

3.740498 

7.214484 

22.363188 

44.241858 

3.740475 

7.214484 

22.363103 

44.241828 

3.740470 

7.214484 

22.363090 

44.241822 

H 

5 

.. 

109.62655 

61.954879 

61.954157 

61.953988 


TABLE II. - COMPARISON OF FREQUENCY PARAMETER RATIOS, (f - P/lj). OF ROTATING BLADES 


FROM BEAM THEORY AND MSC NASTRAN 





r “Pc * 

<p - > ■ 0 , 

(d/b) 

0.05. 




a 

W 1 

Mode 

(a) 

MSC 

NASTRAN 

Beam 

theory 

(non- 

linear) 

Percent 

frequency 

difference 11 

a 

«1 

Mode 

(«) 

MSC 

NASTRAN 

Beam 

theory 

(non- 

linear) 

Percent 

frequency 

difference 6 

0.5 

FI 

F2 

F3 

3.5463 

22.1496 

61.9732 

3.5213 

22.0390 

61.7017 

0.705 

0.499 

0.438 

1.00 

FI 

F2 

F3 

T1 

5.2142 

23.8828 

63.7103 

68.5803 

5.1917 

23.7831 

63.4591 

0.432 

0.418 

0.394 


T1 

SI 

F4 

F5 

68.3947 

69.9391 

121.4804 

201.0141 

70.3203 

120.9065 

199.8640 

-0.545 

0.472 

0.572 


SI 

F4 

F5 

T2 

69.9556 
123.2737 
202.8394 
206 5022 

70.3359 

122.7303 

201.7327 

-0.544 

0.441 

0.546 


T2 

206.2487 

— 






1 0.50 

FI 

F2 

F 3 

n 

si 

F4 

F5 

4.0287 

22.5920 

62.4092 

68.4408 

69.9432 

121.9282 

201.4689 

4.0049 

22.4843 

62.1428 

0.591 

0.477 

0.427 

1.20 

FI 

F2 

F3 

5.7917 

24.6085 

64.4601 

68.6620 

5.7694 
1 24.5132 
[ 64.2176 

0.385 

G.387 

0.376 




T! 

70.3428 

123.5262 

202.5524 

-0.543 

0.427 

0.534 


70.3242 

121.3618 

200.3292 

-0.545 

0.465 

0.566 


SI 

F4 

F5 

T2 

69.9629 

124.0553 

203.6387 

206.6139 


T2 

206.3116 

— 






0.80 

FI 

F2 

F3 

4.6850 

23.2721 

63.0895 

4.6621 

23.1687 

62.8310 

0.489 

0.444 

0.410 

1.50 

FI 
F2 
F3 
1 \ 

6.7182 

25.8909 

65.8156 

68.8120 

6.6958 

25.8031 

65.5890 

0.333 

0.339 

0.344 


T1 

SI 

F4 

F5 

T2 

68.5134 

69.9497 

122.6299 

202.1829 

206.4108 

70.3303 

122.0756 

201.0600 

-0.544 
0.452 
0.555 
- — 


SI 

F4 

F5 

T2 

69.9763 

125.4796 

205.1007 

206.8193 

70.3555 

124.9769 

204.0542 

-0.542 
0.401 
0.510 
1 


*F1 F2 F5 are frequencies In flatwise direction; SI Is first edgewise frequency end 
T1 jnd T2 ere the lowest two torsional node frequencies respect ve y. 

^Percent frequency difference • (fNASTRAN - f Nonl Inear )* lw) ' r NASTRAN- 







TABLE III. - COMPARISON OF LINEAR ANO PERTURBATION FREQUENCIES OF 
POFcnNFD BLADE AT VARIOUS ROTATIONAL SPEEDS 
[B Pc - 15* , f • R - 0, Thickness ratio - 0.05.] 


Q r 

•1 


Mode 


0.3 


0.5 


0.8 


1.0 


2.0 


Beam Theary Results 
(Coriolis Effects Included) 


Linear 

NorI I near 

Percent 

frequency 

change 3 

3.6776 

3.6804 

0.076 

22.1849 

22.1860 

0.005 

61.8467 

61.8475 

0.001 

3.947b 

3.9652 

0.444 

22.4498 

22.4572 

0.033 

62.1113 

62.1169 

0.009 

4.5356 1 

4.6094 

i 1.601 

23.0828 

23.1382 

0.153 

62.7514 

62.7763 

0.397 

5.0138 

5.1444 

j 2.539 

23.6522 

23.7203 

0.287 

63.3360 

63.3821 

0.073 

7.8855 

8.3777 

5.875 

27.9477 

28.3037 

0.444 

67.9947 

68.0724 

0.114 


MSC NASI RAN 
Coriolis effects 
neglected 


fc f Nastran " f Nonlinear! 
I . f Nastran J 


3.7062 

22.2934 

62.1044 

3.9950 

22.5585 

62.3354 

4.6451 

23.2026 

62.8837 

5.1781 

23.7926 

63.4078 

8.3385 

28.3445 

68.0883 


0.696 

0.482 

0.414 

0.746 

0.449 

0.351 

0.769 

0.364 

0,171 

0.&51 

0.304 

0.041 

-0.470 

0.144 

0.023 


^Percent frequency change 


(^nonlinear - f l inear)xlOO/ f nonlinear. 


TABLE IV. 


. FURTHER COMPARISON OF FREQUENCY PARAMETER RATIOS, (f 

. r nftTATI Ahftl SPFFDS 


p/x.J. FROM 


[d/b 


B Pc 

deg 


.05 


0.05 


|0 .05 


0.25 


30 


4b 


45 


45 


Mode MSC 
NASTRAN 


0.9 


0.8 


1.05 


0.8 


1 

2 

3 

1 

2 

3 

1 

2 

3 

1 

2 

3 


4.7719 

23.2410 

62.4392 

4.2769 

22.6971 

61.6110 

4.9615 

23.2508 

61.6563 

4.2598 

12.8590 

22.6460 


Beam 

theory 


4.6937 

23.2313 

62.8783 

4.1026 

22.7218 

62.3800 

4.7622 

23.3468 

62.9538 

4.0236 
b4. 1058 
14.2623 
*>13.9954 
22.7315 
*>22.7209 


Percent 1 

d/b 

“Pc 1 


Mode 

frequency 
difference 3 | 


deg 

W 1 


* U639 

0.042 

0.05 

45 

0.445 

1 

2 

•> 

-0.703 





4.075 



1.0 

1 

-0.109 

0.05 

45 

? 

.1.248 





4.017 

-0.413 

0.05 

60 

0.5 

1 

2 

-2.104 




J 

5.545 




1 

1*3.615 

-’,0.913 

0.25 

45 

1.0 

2 

**-8.837 





-0.378 

*>-0.331 




J 


MSC 

NASTRAN 


J.6511 
22.2438 
61.8224 

4.8337 

23.1064 

61.4499 

3.4708 

',22.2377 

|62.0111 

4.8122 

12.3085 

23.0257 


Beam 

theory 


Percent 
I frequency 
difference 3 


3.6066 

22.2087 

61.8820 

4.6145 

23.2047 

62.8250 

3.3944 

22.1264 

61.8150 

. 4.5128 

*>4.6185 
14.2650 
I bl3.9808 
23.2460 
I b23.2040 


1.219 

0.158 

-0.096 

4.535 

-0.425 

-2.238 

2.201 

0.501 

0.316 

6.222 

*>4.025 

-15.896 

*>-13.587 

-0.957 

*>-0.774 


*P*rC*"t frWOWKy > I Vs,.." ■ Mic N.STftA. 

....... ... !" »<*• “ *" 

calculations) f-Y-0» R-0, 
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TABLE V. - COMPARISON OF LINEAR AND PURTURBATION FREQUENCIES OF ROTATING 
BLADES AT VARIOUS PRECONE ANGLES (CORIOLIS EFFECTS INCLUDED): 


ft/uj . 1.0, « »o . ® " 


"Pc 

Mode 

Low thickness ratio, d/b * 0.05 f 

Hiqh thickness ratio, d/b * 0.25 


Linear 

Non i inear 

Percent 

frequency 

variation 3 

Linear 

Nonlinear j 

Percent 

frequency 

variation 3 

10 

1 

2 

3 

5.1124 

23.7242 

63.4037 

5.1713 

23.7553 

63.4246 

1.130 

0.131 

0.033 

5.0912 

14.1663 

23.7219 

5.1680 

14.1229 

23.7606 

1.486 

-0.307 

0.163 

20 

1 

2 

3 

4.8770 

23.5540 

63.2440 

5.1043 

23.6709 

63.3234 

4.453 

0.494 

0.125 

4.8009 

14.2935 

23.5450 

5.0865 
T 1313 
23.6854 

5.615 

-1.148 

0.393 

:o 

1 

2 

3 

4.4930 

23.2908 

62.9985 

4.9749 

23.5282 

63.1612 

9.687 

1.009 

0.258 

4.3492 

14.477 

23.2718 

4.9316 

14.1553 

23.5526 

11.810 

-2.273 

1.192 

40 

1 

2 

3 

3.9725 

22.9639 

62.6960 

4 . 7c 10 
23.^266 
02.9476 

16.562 

1.555 

C.400 

3.7715 

14.6859 

22.9329 

4.6814 

14.2072 

23.3554 

19.437 

-3.369 

1.809 

45 

1 

2 

3 

3.6659 

22.7880 

62.5344 

4.6145 

23.2047 

62.8250 

20.557 

1.796 

0.463 

3.4464 

14.7907 

22.7509 

4.5133 

14.2477 

23.2316 

23.639 

-3.811 

2.069 

50 

1 

2 

3 

3.3315 

22.6109 

62.3724 

4.4371 

23.0697 

62.6942 

24.917 

1.989 

0.513 

3.1023 

14.8916 

22.5677 

4.3129 

14.3002 

23.0910 

28.069 
-4.136 1 

2.266 | 


^Percent frequency variation - (“nonlinear - “1 i near ) “von 1 inear . 



TABLE VI. - EFFECT OF LINEAR AND NONLINEAR CORIOLIS FORCES. ANO VARIOUS NONLINEAR 
TCRNS OH THE FREQUENCY PARANETER RATIOS FOR UNTWISTED THIN BLADES 

(») d/to . 0.05. , . 0*. T - 0\* . 0 


*Pc 

0 

"1 

Node 

Solution of linear equations 

Perturbation solution of nonlinear equations 

L_ 




Coriolis 

Coriolis 

Full 

Linear 

— 

Nonlinear 

A 1jk • ° a 

D,j k - °? 

Df-jit - 0 




forces 

forces 

nonl Inear 

Coriolis 

Coriolis 

Eljk - Ob 

E jk - 0 




included 

ignored 

equation 

forces 

forces 


P jk - o e 







ignored 

ignored 



QijSE - ° c 



1 

3.9476 

3.9480 

3.9652 

3.9654 

3.9650 

3.9593 

3.9535 

3.9535 


o.s 

l 

22.4498 

22.4498 

22.4572 

22.4572 

22.4572 

22.4563 

22.4507 

22.4507 



3 

62.1113 

62.1114 

62.U69 

62.1169 

62.1173 

62.1170 

62.1115 

62.1117 



1 

4.5356 

4.5366 

4.6094 

4.6096 

4.6087 

4.5852 

4.5608 

4.5608 

5* 

0.8 

2 

23.0828 

23.0829 

23.1182 

23.1182 

23.1185 

23.1142 

23.0873 

23.0874 



3 

62.7514 

62.7515 

62.7763 

62.7763 

62.7807 

62.7792 

62.7525 

62.7531 



1 

5.0138 

5.0155 

5.1444 

5,1443 

5.1431 

5.1021 

5.0689 

5.0589 


1.0 

2 

23.6522 

23.6524 

23.7203 

23.7202 

23.7210 

23.7131 

23.6611 

23 6611 



3 

63.3360 

63.3361 

63.3821 

63.3820 

63.3936 

63.3909 

63.3383 

63.3j94 



1 

3.5553 

3.5576 

3.6482 

3.6 98 

3.6476 

3.6180 

3.5872 

3.5672 


0.5 

2 

22.2252 

22.2256 

22.2616 

22.2618 

22.2616 

22.2574 

22.2299 

22.2299 



3 

61.9077 

61.9078 

61.9351 

61.9352 

61.9364 

61.9350 

61 9083 

61.9094 



1 

3.6141 

3.6199 

4.1026 

4.1050 

4.0991 

3.9591 

3.7994 

3.7994 

15* 

0.8 

2 

22.5196 

22.5206 

22.7218 

22.7221 

22.7230 

22.7050 

22.5476 

22.5478 



3 

62.2344 

62.2348 

62.3800 

62.3800 

62.3976 

62.3943 

62.2388 

62.2448 



1 

3.6659 

3.6751 

4. 6145 

4.6164 

4 .6069 

4.3631 

4.0621 

4.0621 


1.0 

2 I 

22.7880 

22.7895 

23.2047 

23.2045 

23.2080 

23.1816 

22 8500 

22.8504 



3 j 

62.5344 

62.5349 

62.8250 

62.8246 

62.8761 

62.8784 

62.5451 

62.5576 



1 

3.0412 

3.04 50 

3.0412 

3.0450 

3.0412 

3.0412 

3.041? 

3.0412 


0.5 

2 

21.9636 

21.9643 

21.9636 

21.9643 

21.9636 

21.9636 

21.9636 

21.9636 



3 

61.6719 

61.6722 

61.6719 

61.6722 

61.6719 

61.6719 

61.6719 

61.6719 



1 

2.1029 

2.1096 

2.1029 

2.1096 

2-1029 

2.1029 

2.1029 

2.1029 

>0* 

0.8 

2 

21.8524 

21.8542 

21.8524 

21.8542 

| 21.8524 

21.8524 

21.8524 

21.8524 



3 

61.6324 

61.6331 

61.6324 

51.6331 

! 61.6324 

61.6324 

61.6324 

61.6321 


1.0 

1 

2 

unstable 

unstable 

unstable 

unstable 

unstable 

unstable 

unstable 

unstable 



3 


— 


1 



1 




^N an) i near terns 
b Non) inear terms 
c Nonl inear terms 


due to foreshortening are iqnored. 

arising from (Tv ') ' and ( Tw 1 ) * are ignored. 

associated with extentional deformation are ignored. 


TABLE Vi. - EFFECT OF LINEAR ANO NONLINEAR CORIOLIS FORCES, AND VARIOUS NONLINEAR 
COUPLING TERMS ON THE FREQUENCY PARAMETER RATIOS FOR PRETWI3TE0 THIN BLADES 

(b) d/b . 0.05, , - 0*. Y - 30\ R - 0 


ape 

Q 

-1 

Node 

Solution of linear equations 

Perturbation solution of nonlinear equations j 

Coriolis 
forces 
i nc 1 uded 

Coriolis 

forces 

ignored 

Full 

nonlinear 

equation 

Linear 

Coriolis 

forces 

ignored 

Nonl inear 
Coriolis 
forces 
ignored 

A 1jk ■ 0 

o o 

a • 

-K 

P'Jk - 0 

§ 1 jk - 0 
p ijk “ 0 

Qijk - o 



1 

3.9540 

3.9551 

3.9716 

3.9772 

3.9711 

3.9657 

3.9600 

3.9599 


0.5 

2 

20.0667 

20.0621 

20.0716 

20.0688 

20.0733 

20.0715 

20.0664 

20.0664 



3 

59.0449 

59.0442 

59.0499 

59.0497 

59.0504 

59.0499 

59.0449 

59.0450 



l 

4,5355 

4.5388 

4.6099 

4.6107 

4.6079 

4.5855 

4.5613 

4.5613 

15 

0.8 

2 

20.6248 

20.6132 

20.6488 

20.6460 

20.6563 

20.6477 

20.6232 

20.6232 



3 

59.6203 

59.6184 

59.644? 

59.6447 

59.6470 

59.6441 

59.6199 

59.6202 



1 

5.009? 

5.0147 

5.1416 

5.1417 

5.1374 

5.0986 

5.U5S7 

5.0557 


1 .0 

2 | 

21.1262 

21.108? 

21.1731 

21.1724 

21.1869 

21.1700 

21.1225 

21.1225 



3 

60.1438 

60.1408 

60.1903 

60.192? 

60.1962 

60.1904 

60.1428 

60.1434 



1 

3>5b^' i 

3.56?? 

3.6519 

3.6578 

3.6499 

3.62)6 

3.5913 

3.5913 


0.5 

2 

19.8932 

1 g 

19.9169 

19.8912 

19.9252 

19.9171 

19.8922 

19.8922 



3 

58.86fl 

vu.ob 94 

58.8890 

58.8857 

58.8915 

58.8889 

58.8647 

58.8653 



1 

3.6066 

3.6265 

4.0952 

4.1046 

4.0841 

3.9503 

3.7928 

3.7928 

45* 

0.8 

2 

20.1907 

20.1035 

20.3266 

20.2884 

20.3700 

20.3268 

20.1820 

20.1821 

J 


3 

59.1652 

59.1506 

59.2997 

59.2983 

59.3148 

59.3035 

59.1620 

59.1654 



1 

3 . 6481 

3.6791 

4.6005 

4.6089 

4.5765 

4.3454 

4.0477 

4.0477 


1.0 

2 

20.4614 

20.3265 

20.7467 

20.7138 

20.8328 

20.7483 

20.4387 

20.4388 

— 

— 

3 

59.4402 

59.4172 

59.7200 

59.7250 

59.7531 

59.7361 

59.4320 

59.4388 


Mode Solution of linear equations 


Coriolis Coriolis 
forces forces 
included ignored 


3.9390 

14.1019 

22.4485 
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TABLE VII. - EFFECT OF LINEAR AND NONLINEAR CORIOLIS FORCES, AND VARIOUS NONLINEAR 

rfinoi T nr Trnur nu tiic rnr^itrury n«n«urtpft »»t, T i>f 

l*i»u iuImu wii m*L 1 r mhl i l rv nniiLU rur\ irucfk OCMUCJ 


(a) d/b - 0.25, , . 0", y - O'.ti - 0 


Perturbation solution of nonlinear equations 


Full Linear 
nonlinear Coriolis 
equation forces 


Nonlinear 
Coriol i s 
forces 
ignored 


69 


14.5808 

22.6865 


14.8973 

23.1473 


14.0964 

22.4550 

14.0947 

22.4494 

14.0948 

22.4497 

4.5707 

14.1351 

23.1111 

4.5464 

14.1268 

23.0842 

4.5458 

14.1270 

23.0852 

5.0816 

14.1654 

23.7090 

5.0388 

14.1486 

23.6570 

5.0376 

14.1488 

23.6579 

3.5690 

14.2234 

22.2469 

3.5387 

14.2150 

22.2194 

3.5387 

14.2155 

22.2220 

3.85^2 

14.3878 

22.6747 

3.6990 

14.3300 

22.5167 

3.6989 

14.3331 

22.5308 


Solution of linear equations 


Perturbation solution of nonlinear equations 


Coriolis 

forces 

Included 

Coriolis 

forces 

ignored 

Full 

nonlinear 

equation 

Linear 

Coriolis 

forces 

ignored 

Nonlinear 

Coriolis 

forces 

ignored 

A ijk “ 0 

F ijk ‘ 0 
E 1jk - 0 

Ojjk - 0 

E ijk „ 
j!l jk * 0 
Qijk * 0 

3.9439 

3.9538 

3.9649 

3.9705 

3.9611 

3.9576 

3.9518 

3.9517 

13.2169 

13.1901 

13.2064 

13.1937 

13.2191 

13.2120 

13.2098 

13.2099 

24.0327 

24.0341 

24.0394 

24.0408 

24.0388 

24.0381 

24.0331 

24.0335 

4.5074 

4.5369 

4.5995 

4.6048 

4.5800 

4.5687 

4.6446 

4.5441 

13.3419 

13.2712 

13.2948 

13.2944 

13.3524 

13.3169 

13.3061 

13.3062 

24.5895 

24.5933 

24.6244 

24.6279 

24.6198 

24.6164 

24.5923 

24.5936 

4,9606 

5.0122 

5.1294 

5.1273 

5.0888 

5.0750 

5.0325 

5.0313 

13.4530 

13.3393 

13.3642 

13.3926 

13.4734 

13.4016 

13.3805 

13.3804 

25.0954 

25.1018 

25.1669 

25.1711 

25 1554 

25.1492 

25.1022 

25.1037 

3.5016 

3.5658 

3.6075 

3.6553 

3.5907 

3.5708 

3.5410 

3.5410 

13.3418 

13.1097 

13.2916 

13.1247 

13.3530 

13.3197 

13.3084 

13.3089 

23.8329 

23.8359 

23.8630 

23.8664 

23.8616 

23.8578 

23.8337 

23.8367 

3.4631 

3.6246 

4.0136 

4.0863 

3.9205 

3.8410 

3.6897 

3.6896 

13.6592 

13.0703 

13.4067 

13.1819 

13.7269 

13.5344 

13.4611 

13,4639 

24.0827 

24.0918 

24.2659 

24.2760 

24.2461 

24.2301 

24.0899 

24.1062 

3.4281 

3.6766 

4.4952 

4.5530 

4.2956 

4.1909 

3.9127 

3.9126 

13.9447 

13.0318 

13.4582 

13.3039 

14.0933 i 

13,6823 

13.5075 

13.5130 

24.3105 

24.3268 

24.7071 

24.7230 

24.6541 

24.6313 

24.3319 

24.3640 






























TABLE VIII. - CrrCCT or LINEAR AND NONLINEAR CORIOLIS FORCES, AND VARIOUS 
NONLINEAR TERMS ON THE FREQUENCY PARAMETER RATIOS 


[Presented as percent frequency variation based upon full nonlinear 
equation solution, (f n onl Inear - f reference) xl0 °/f nonlinear-] 


(j) Bp c * 15 # , f • y • 0*, R ■ 0 


a 

Mode 

All nonlinear 

Linear Coriolis 

Nonlinear Coriolis 

A 1 Jk * 0 

D 1jk * 0 

-1 


terms Ignored 

forces Ignored 

forces ignored 

Eljk - 0 

Thin blade: (d/b) - 0.05 

0.5 

1 

0.4439 

-0.0050 

0.0050 

0.1488 

0.2951 

0.8 

1 

1.6011 

-0.0043 

0.0152 

0.5250 

1.0544 

1.0 

1 

2.5387 

0.0019 

0.0253 

0.8223 

1.6620 

0.5 

2 

0.0330 

0.0 

0.0 

0.0040 

0.0289 

0.8 

2 

0.1531 

0.0 

-0.0013 

0.0173 

0.1337 

1.0 

9 

0.2871 

0.0004 

-0.0030 

0.0304 

0.2496 

0.5 

3 

0.0090 

0.0 

-0.0CD6 

-0.0001 

0.0087 

0.8 

3 

0.0397 

0.0 

-0.0070 

-0.0046 

0.0379 

1.0 

3 

0.0727 

0.0002 

-0.0181 

-0.0139 

0.0691 

j Thick blade: (d/o) - 0.25 

0.5 


0.5303 

-0.1263 

0.0884 

0.1843 


HO 

1 

1.9865 

-0.0956 

0.3S12 

0.6607 

1.1889 

1.0 

i 

3.2463 

0.0526 

0.7283 

1.0419 

1.S753 

0.5 

2 

-0.0809 

0.0795 

-0.0916 

-0.0419 

-0.0298 

0.8 

2 

-0.3621 

-0.0262 

-0.4146 

-0.1658 

-0.1070 


2 

-0.6739 

-0.2407 

-0.7787 

-0.2803 

-0.1614 


3 

0.0370 

-0.0031 

0.0053 

0.0080 

0.0330 

0.8 

3 

0.1817 

-0.0035 

0.0316 

0.0450 

0.1613 

1.0 

3 

0.3485 

-0.0021 

0.0636 

0.0872 

0.3064 


(b) Bp c =45 ,9 m y m 0*, R « 0.] 
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TABLE IX. - EFFECT OF LINEAR AND NONLINEAR CORIOLIS FORCES, AND VARIOUS 
NONLINEAR TERMS ON THt FKEgUENCV PARAMETER RATIOS 

[Presented as percent frequency variation based upon full nonllncir 
equation solution.]) 

(a) »p c - 15*, v - 0*, R • 0, y - 30 

[ Mode | All non'llnearl Linear Coriolis) Monl1ne 4 r Corfolls| A 1jk - « Kjk 
terms Ignored forces Ignored force s Ignoreo 

~ Thin blade: (d/b) . 0.05 

n 1 o~ 4432 I -0.0151 I 0.0126 0.1486 I 0.2921 

1 ?* 5751 oC US US 


0.0003 

-0.0008 

-0.0032 

Thick blade: (d/b) 

-0.1412 
-0.1152 
j 0.0409 


-0.0058 

-0.0142 

-0.0207 


(b) Bp c - 45 , v - 0 , R • 0, y - 31 

Mode I All nonlinear I Linear Corlollsl Nonlinear Coriolis A 1jk 
terms ignored forces Ignored forces ignored 


0.0126 

0.1486 

0.2921 

0.0434 

0.5293 

1.0543 

0.0817 

0.8363 

1.6709 

-0.0085 

0.0005 

0.0259 

-0.0363 

0.0053 

0.1240 

-0.0652 

0.0146 

0.2390 

-0.0009 

0.0 

O.OOOC 

-0.0047 

0.0002 

0.0407 

-0.0098 

-0.0002 

0.0789 

0.25 



0.0958 

0.1841 

0.3304 

0.4240 

0.6696 

1.1936 

0.7915 

1.0606 

1.8891 

-0.0962 

-0.0424 

-0.0258 

-0.4333 

-0.1662 

-0.0850 

>0.8171 

-0.2799 

-0.1220 

0.0025 

0.0054 

0.0262 

0.0187 

0.0325 

0.1304 

0.0417 

0.0664 

0.2531 


2.5329 

11.9310 

20.7021 


Thin blade: (d/b) 

r -0.1616 

-0.2295 

-0.1826 


0.0056 

0.0024 

-0.0084 

Thick blade: (d/b ) 


-0.0417 

-0.2135 

-0.4150 

-0.0043 

-0.0255 

-0.0554 



1.6594 

7.3843 

12.0161 


0.0002 
-0.0064 
j -0.0270 


2.9356 

13.7159 

23.7387 


1.8434 

8.0701 

12.9583 

-0.1264 
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Figure 3. - Effect of pretwist, precone and setting angle variations on the fun- 
damental mode frequency parameter of rotating, thin, blade: Perturbation 
solution: d/b ■ 0.05, Q/u^ * 0.8, L/b * 10, L ■ 5.0 in. 
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